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On Harmonic Functions.* 

By R. A. Harris. 



§1. 
This paper constitutes an elementary study in harmonic functions, V, and 
what might be called quasi-harmonic functions, V, where 

Oi 2 17* C\ 2 T7" 

3^- + ^yr+""=°» or V 3 F p = 0, (1) 

and 

3 2 V /3 8 V \ 

a^-CaF"+----; = ' or v2V ^ = ' < 2 ) 

jp denoting the number of independent variables. While its object is to show 
that many solutions of (1) and (2) can be found by elementary methods and to 
point out some of the connections which exist between various solutions, the 
fact that many, and perhaps all, of the solutions may, in a sense, be regarded 
as an outgrowth or extension of those harmonic functions, X, Y, with which the 
theory of functions has to deal, will not be lost sight of. 

The partial differential equations (1) and (2) will not be transformed into 
their equivalents in spherical or cylindrical coordinates even when solutions 
involving Legendre's coefficients, BessePs functions, etc., are contemplated. In 
adapting equation (2) or any of its solutions to physical problems, the quantities 
x, y, z are generally to be replaced by xt, x, and y, respectively. 

§2. 
The lemmas or propositions given below are, with few exceptions, self- 
evident or nearly so. But inasmuch as their wording may not always be exact 
and such as to cause their meaning to be readily apprehended, brief demon- 
strations and qualifying statements are added. 



* An abstract of the principal parts of this paper was read before Section A of the American Association 
for the Advancement of Science, 1898. See BvXl. Am. Math. Soc, Vol. V (1898), pp. 96-98. 
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(i) If a solution of any linear partial differential equation having constant 
coefficients, but not a term of order zero, be differentiated (or integrated) with 
respect to one or several of its variables, the result is a solution of the given 
equation. 

For, when several operators like ^— , » — ^— , .... are applied to any 

*C 3C y 

quantity or solution the order of application is immaterial ; the result of any 
number of such operations on the given differential equation is, of course, a 
new differential equation. So if a partial derivative of the function constituting 
the given solution be made, by transposing the operators, to occupy the same 
position in the new differential equation as did the given solution in the original 
equation, the result must still be a true equation. 

Similarly for partial integration ; but it is understood that in this case the 
number of terms in the power-series representing the given solution developed 
in powers of the variable in question, is generally infinite. 

(ii) If in any solution of an equation of the kind just referred to, any two 
of the variables be replaced by their equivalents in circular coordinates, the 
result of differentiating (or integrating) the solution with respect to either circular 
coordinate is a solution of the given equation. 



For, suppose rj = y + s/ — lz, £ = y — \/ — lz, and so 

d _ d dy d dz 1 / d . d \ 

dri dydri dzdvi ~*\dy "~ dz)' 

3 _ x ( A _l / t d \ 

dl;- 2 \dy + v - l dz')' 

and so, by lemma (i), these indicated differentiations yield solutions. Similarly 
for coordinates such that, say, y + z = ri,, y — % = %/, which maybe called 
hyperbolic coordinates. 

The same is evidently true if any other linear substitutions be made for 
the variables. 

(iii) If in any solution of (1) the letters x, y, z, be interchanged at 

will, the result is a solution. The same is true for all variables after the first 
in a solution of (2). 

For, the operators in (1) are symmetrical in all variables, and in (2), in all 
but the first variable. 
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(iv) If in any solution, V p , of (1), y, z, . . . . be replaced by ± V — 1 y, 

± */ — 1 z, . . . . , the result is a solution, V p , of (2) ; and conversely. 

For, equation (1) must still be a true equation provided the substitution 

be made in the operators as well as in the original function Y v . But the 

32 32 32 32 

operators ~— 5 , ~— 3, • • • • become — ^ — 5 , — «— 5, ..... and so the transformed 

r oy 2 az z dy z oz 2 ' ' 

function must satisfy (2). 

(v) Given any solution of (1) or (2) in the form 

/o(y,a,----)x(*)+/.(y,a,....)x"(*)+/*(y,«,---Oar (»)+■■••, (3) 

and such that the entire coefficients of %"(%), % IV («), •••• in (1) or (2) are 
severally zero, then is 

My,z,....)Mx)+U{y,z,....)W{x)+f i < < y,z,....)V*{x) + .... (4) 

a solution, -^ being an arbitrary function. 

This means that the given solution when substituted in (1) or (2) satisfies 
the same irrespective of the nature of x> an & so X ma y De replaced by $. 
To say that the coefficient of % (v) in (1) or (2) shall be zero, is equivalent 
to saying that 

/,-2 (y, a, • • • •) ± (,ap + a^+ • • • •)/*(&, z > ••••) = °- ( 5 ) 

Of course this is likewise true for the alternate terms; that is, for the terms 
having odd subscripts and indices. 

(vi) If in any solution of (1) or (2) referred to in lemma (v), x, y, . . . . 
be replaced by x ± x', y ± y', . . . . , where a/, y', ■ ■ • • are constants, the result 
is also a solution of the given differential equation. 

For, let x, = x d= x', y,=zy ±y' ; the solution will be the same 

function of x, , y, , . . . . as it originally was of x, y, The operators 

~— . =s— , • • • - may be written ^ — , « — , But expression (5) being 

dx' oy ox, oy, r v 

identically zero, it will remain so when x, , y, , .... replace x, y, .... 

throughout. 

(vii) If for the variables in equation (1) we make an orthogonal sub- 
stitution, the result is still a true equation. 

For, in V, replacing x, y, . . . . by their equivalents in a second set of 
50 
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variables x,, y,, . . . ., the form of V will generally be altered, but its value 
will, of course, remain the same as before. Consequently the operator 

d 2 . d 2 



dx 2 dy 



+ ■ 



But the 



applied to this new expression will still give a result equal to zero. 

effect of the operator 

d 2 d z 

dx~f + dy~f + "" 

upon either form of V is not zero unless the two operators are equal or propor- 
tional. The first operator when expressed in terms of x,, y lf . . . • by aid of the 
matrix of transformation shows that for any orthogonal substitution the two 
operators are equivalent, ignoring a constant factor, the modulus of the sub- 
stitution. For suppose 



x =za 1 x, + b 1 y,+ . 
y = a z x, + b z y,+ . 



and so - 



then 



dx 



d dx/ 
dx t dx 



9 dy, 
dy, dx 



+ ^^r + 



x, = a 1 x + a z y + 

y, = bin + hy + 



1 ox, oy. 



■ > 



and finally, 

d 2 , d 2 



dx 2 
wherein 



+ gp +...-= (af + a 2 + 



■)^ + (hl + bl+....)^ + 



«! + «!+• 



5S + 6|+....= 



Similarly, if only a part or group of the variables of (1) be transformed 
into a group of corresponding variables by means of a uni-modular orthogonal 
substitution, the equation will still hold true. This remark is evidently appli- 
cable to the second part of equation (2). 

(viii) Equation (1) is satisfied by making V equal to 

d p - 2 log p ... 1 



or 



,p-3> 



dp p - a p J 

where p denotes the number of variables in V and where 

p = V x 2 + y 2 + 



(«) 
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For, substituting l/p p ~ z in (1) we have the sum of p quantities like 

(p-2){-p-»+pz*p-*-'), 
except that x z is replaced by y z , z z , . . . . in succession; and this sum is zero. 
Similarly, or by lemma (iii), 

&=» ( 7 ) 



where a = */ x z — (y 2 + ■•■•)> 

is a solution of equation (2). 

(ix) If H n denote an homogeneous solution of (1) of degree n, then is 



E n 



-Zn + p — Z 



(8) 



likewise a solution, n being general. 

d z 
For, operating upon H n /p Zn+p ~ z with =— g, we have 

_J 8»J5k 2g(2n+f> — 2) 8.ff w 2n+ j>+ 2 „• 

pZn+p-Z Q x » p Zn + p ~daT p Zn + P n 

x z (2n + p — 2) (2n + p) ff 
+ zn+p+z ■**»• 

r 

3 2 
Replacing a; by y, we have the result of the operation ^— 3 ; and so on. But H n 

d s H 
being a solution, the result of adding together the terms involving „ 8 " , 

O fly 

-^ — ^, .... must be zero. Because H n is homogeneous and of degree n, the 

sum of the second terms is * — aJ+p ^«- The sum of the third terms 

r 

is — — — zwtp &n 5 the sum of the fourth terms, since p z = x z + y z + • • • • , 

P 

is — 8w+p — #„ • Hence the total is zero. 

r 

Conversely, if (8) be a solution and H n an homogeneous function, Z7 W must 
be a solution. 

Similarly, or by lemma (iii), if H n denote an homogeneous solution of (2), 

^Zn+p—Z > \"/ 

where 

a = V x z — (y z + z z •+- ), 

is a solution. 
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§ 3. Application of the Foregoing Lemmas to Spherical Harmonic Functions. 

If the solutions (6), (7) be differentiated n times with respect to x, we 
obtain solutions of the forms 



where i^T — J, J^f— ) denote algebraic integral functions of the «-th degree 

in aj/p, aj/c. 

The truth of this statement becomes evident upon performing a few differen- 
tiations and referring to lemma (i). The converse of (ix) shows that the 
homogeneous functions 

fj-.(|), «•*.(*) (id 

are also solutions of equations (1) and (2). These solutions evidently have 
circular symmetry about the as-axis; they are termed zonal harmonics of the 
n-th order while (10) are zonal harmonics of order — n — p + 2. 

If the solution l/p p-iJ be differentiated with respect to a circular coordinate, 
£, we have 

P P ~ 2 = / 1)n (P- *)P (P + 2) . ■ ■ ■ Q + n - 1) n n h2 x 

where yi and £ are conjugate imaginaries. 

Similarly, if the solution of (2), l/cr p-2 , be thus differentiated, we have 

3» 1 

T p - 2 - ( , 11 n (p-%(P+2)....( j) + n-l) q" n ,x 

\~T~ l ) ™ ~v + 2(n-l) • v 1 ^ 



3/» v ' > 2 n a P + 2(n-l) 

If >7 /; £, denote hyperbolic coordinates such that 

vi, = a; + y (or s or «?, . . . . ), 
%t = x — y (or aorw, ••••), 
the n-th derivative with respect to £, is 

2» * _ 

/ 1 x ra (i>- 2 )P(i>+ 2 )---(jP + ^-l) >?? r , 4 x 

V ' on -p + 8(n-l) ' V / 



OfP-8 



a?r 
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In all of the above cases 37 and £ (or yj, and £,) can be interchanged at will. 
By (ii) and (vii), the expressions involving p are solutions of (1), and those 
involving a, solutions of (2). 

If the solutions 

and —r— o , or 



P ff (£C 2 rfc C?) nr 

where <?, — y % -\- z % + ..... be integrated one or more times with respect to x, 
we obtain solutions of (1) and (2). 

Making p = 3 and performing one integration, we obtain from the solutions 

rfc - , ± - , the solutions 

1 i— r + g l lnff «r+ai 

The former may be written 



log 7 ^, ±log°±*. (15) 

2 & r — a;' 2 & <r — 05 v/ 



| lo ST^= r °«o(^ (16) 

where ^ = x/r. Another integration gives the solutions 

1 , r + x 1 , c + x /ir ,\ 

The former may be written 

Solutions like these are zonal harmonics of the second kind. 

Zonal Harmonics, p = 3. By lemma (viii) or by direct substitution, we 
see that 

v 1 1 l y z + z z 1.3 (y» + s»)« 1.3.5 (y' + s 8 )' , hq , 

Vs ~Vx 2 -\-y^ + z 2 ~ x 2 x 3 i "2.4 a 6 2.4.6 x" t- • ■ • • V lw ; 

satisfies equation (1). Moreover it is a solution of the form referred to in (v), 
%{x) being 1/x. Therefore, by (vi), 

*/ (x — x') 2 + (y — y'f + (a - a') 2 ^ 

satisfies equation (1), as can be directly verified. The first radical denotes the 
distance of a variable point from the origin ; the second, the distance of this 
point from a fixed point whose coordinates are x', y', zf, both points being 
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referred to the same coordinate axes. Suppose the fixed point to lie upon the 
x- (or a;'-) axis; i. e., suppose y' = z' = 0. The distance of x, y, z from 
x', 0, is V (x — x'f + y z -\-z z or Vr 2 - 2xr l + r' z , where r z = x z + y z + z z , 
r' z = x' % + 2 + 2 = x lz . Expanding by the binomial theorem and arranging 
the result in powers of r/r' or r'/r, we have 

Vr z -2xr' + r' z r ,L_ 2 ^ ^ ,Jl-2^ + ^ ' ' 

= ^[P + P 1 ^ + P a ^ a +....] (22) 

1 T r' r' 2 1 

= yL p ° +Pi r + P2 ^ +••••]• ( 23 > 

In these developments it is evident that P denotes a rational integral 
algebraic function of x/r. But (22) or (23) satisfies Laplace's equation (1), 
regardless of the value assigned to the constant r' or to the ratio r'/r. Conse- 
quently since (22) and (23) are arranged in powers of this arbitrary magnitude, 
each of their terms is likewise a solution. In other words, 

r»P,(i) OT r— '?,(£), (24) 

where n = 0, 1, 2, . . . . , is a solution. [Cf. lemma (ix).] 
We may write r n P n f~j in the form 

A n x n + B n x n ~ z r z + C n x n - i r i -\- ...., (25) 

or K n<0 x n + K ny% x n ~ z r z + K nA x n -* r" + ..... (26) 

or 2 K n ,x n ~'r'. (27) 

y = 0.2,4 

Expressions for A n , B n> C n , . . . . may be obtained from the expansion of (21) 
in powers of r/r 1 , by picking out the entire coefficient of the «-th power. The 
sum of the coefficients entering this entire coefficient may be shown to be unity. 
Laplace's Coefficients. Returning to solution (20), in which x', y', z' may 
refer to any fixed point in space, we have, since r' z = x' z -f y' z + z' z , 

1 



V (x — x') z + (y — y'f + (z — z') z 
1 1 



, / ^ xx' + yy' + zz' r r z I ^ xx' + yy' + zz' r' , r' z ' 

\ rr' r'^ r' z \ r~r' r r> 



(28) 
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These expanded by the binomial theorem and arranged in powers of rjr' 
or r'/r may be written 

^[ y o+3 7 i^ + F 2 ^+....], (29) 

f[r.+ r 1 £ + r,£+....]. (so) 

It is evident that Y denotes the same rational integral function of 
xx' + yy' + zz' that p (L egendre > 8 coefficient) denotes of x/r. The Y's are 

Laplace's coefficients. The expression ^ is the cosine of the angle 

at the origin between the point x, y, z and the point x 1 , y', z'. Since the value 
of r'/r is arbitrary, each term of the above expressions is a solution of (1). 

TesseraZ Harmonics, Associated Forms. Any unimodular orthogonal sub- 
stitution does not alter the value of r, and the result is by lemma (vii) a solution 
of equation (1). 

In the solutions r n P n and r~ n ~ 1 P n , r and so r n and r~ n ~ l are unaltered 
in value if we put x = l t x, + m x y, + «i z, , where l 1} m lf n x are quantities 
subject to no other condition than 7f + m\ + n\ = 1. Therefore, if l x = 1, 
then m 2 + n\ = or % = \/ — 1 wij . Let mj = c^ ; then % = V — 1 c^ and 
as = a?, + a a {y, + V — 1 z,). In substituting this value for x in expression (25), 
we have to consider only the explicit x since r remains unchanged. The result 
of the substitution arranged in powers of a 1 (y, + V — 1 z,) will be a series of 
expressions of the form 

<*! (y + *f — i 2 ) s [» ( w — 1)- ••■(»— s + l)A*"~ s 

+ (n — 2)(« — 3).... (n — s— 1) -B„ a"-"- 2 r 2 - (31) 

+ (w — 4)(n — 5).... {n — s—S)O n x n - s - i r i +....], 

wherein all subscripts have been omitted from the variables. Bach expression 

of the above form is a solution of the equation (l) ; for, the entire expression (25) 

is still a solution after the above transformation has been made, whatever be the 

value of «!• Consequently, a t being arbitrary, the entire coefficient of any 

power of <*! must be a solution. Further, the real part will be one solution, 

the imaginary part another. 

By writing x = r cos 0, =rp, y = r sin cos $, z = r sin 6 sin <£>, expression 

(31) when divided by af becomes 

d* p 

r n sin 8 (cos sty -f V — 1 sins 4>)-j— r> ( 32 ) 

which may be written 

r n (cos s $ + V — 1 sin s #) T™. (33) 
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Expression (27) may have its y z + z z replaced by the equivalent, vi£, in 
circular coordinates and, by lemma (ii), the result of one or more differentiations 
with respect to £ will be a solution. 

A 2 K„x n -'r" = v J) \K^x-*r*-\ (34) 

' x = 0,2,4 ^ = 2,4 

By making s = 1 in (31), the latter becomes, to a factor, 

y, 2 {n — v)K n ^x n - l -"r\. (35) 

■- = 0,2,4 

Upon replacing n by n — 1 , v by v — 2, this becomes 

n 2 (n — r + l) £"»-,,,_„ as n -r-». (36) 

•• = 2,4 

Here are two solutions (34), (36) of equation (1), both of which are identical 
so far as the variables are concerned. Consequently, the coefficients of like 
powers of the variables must have a constant ratio to each other. 

• _ H— — - -g j^ - = constant, = i. (37) 

• ' 2(n— v + 1) ^_i,,_, '2' v ' 

as any special case will show. This furnishes a ready means of deriving P n 
from P„_i ; for 

-S-n, v = -O-n-l, x-2 > (38) 

= - ( "- 1) 7 ( "- 2) ^,-,-,. (39) 

Example. Given P x = — or /it, to determine P 2 , P s , P it Here 



■Ki.o = i; 



r>- 2 — 2 + 1 J-?- 1 

- • -"-2, 2 g *' ° — 2 ' 



The sum of all coefficients in a P function is to be taken equal to unity. 

- — __i 

2 r 2 ~ 2 ' 

Similarly, 



. p — 

• -^2 — o ^2 



5 A3 3 3 X, 

8 r 4 — 8 r 8 + 8 



P — - _ _ 

•*3 — « ~» ~ 



Other matters concerning spherical harmonics will be found in §§ 4 and 6. 
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§ 4. Applications of Lemma (v). 

The solutions l/p p ~ 2 , lja p ~ s of equations (1) and (2) become when expanded 

in powers of x and c,, where cf = y % + a 2 + , 

f ,_ 1V ( j> - 2)j> ( p + 2) . . . . (p + 2 y - 4) 1 c," f . 

In particular, if ^> = 3, 

l_l_ly 2 + 2 2 1.3 (y 2 +/) 2 _ 1.3.5(y 2 + s 2 ) 3 , 

r~ x 2 a; 3 "'"2.4 sb* 2.4.6 x 7 ^ ' v ' 

\i % (x)=\, then % »(x) = |-, z lv (*)=^, •■••! 

/. ^ = ^^)-^ ! (2/ 2 + « 2 )^(^) + ^^(^ + ^)^ IV W- (42) 

Let ^ (^) = sm h * > t nen 4'" (») = sinh x ; 4' 1V (») = sinh a;, ; 

.-. Binh»[l-.J,(y" + B")+ j^i (y' + «')'-••••], (43) 



or 



sinh x Jq{-ob), (44) 

where w 2 = y z + z 2 and J denotes Bessel's function of order zero, is a solution 
of Laplace's equation (1). 

If 4 (x) = x n , then ty (as) = m (n — 1) ff M - 2 , 

^ IV (a;) = n (n — 1) (« — 2) (w— 3) a;"- 4 , ; 



• ^H ._■ "(*-!) y 8 + a 8 I n(n-l)(»-2)(n-3)(y 2 + g a ) a 1 , , 

• • * L 2 1 a; 2 "*" _ 2 2 4 Z x 4 * ' * ' J K } 

is also a solution. Putting x = rjx, this may be written 
or 



^„n[\ wfri-1) l-/* 2 ■ n(n-l)(n-2)(n-3)(l-^) 2 _ 1 



since, by writing r 3 — a; 2 for y % + g 3 , (45) is of the form (25) and r n P n (fi) 
satisfies equation (1). 

[It may be well to note that the general expression (45) for a solid zonal 

harmonic can be obtained by stai'ting with the solution a: 2 — - (y 2 + z 2 ), where 

% (x) = a; 2 , proceeding to one of the third degree by making 4 (#) = * 3 , and 
then to a solution of the fourth degree by making 4* («) = »*, an d so on. 
51 
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The coefficient of the added term in the solution of the fourth degree, viz., 
Ki(y z + z 2 ) 2 ^ {x), is found by substituting the solution in equation (1). In 
general, it will be found, upon equating to zero the entire coefficient of x n ~ r y"~ 2 
in (1), that v 2 K v = — K v _ 2 . Hence, any solid zonal harmonic may be written 

y 2 + z 2 d 2 (x n ) (y 2 + z 2 ) 2 g*(g") 






i 



2 2 dx 2 T 2 2 .4 2 dx 
which agrees with (45).] 

Replacing n by — n — 1, we obtain the solution 
1 r (n + !)(* + 2) l-(i 2 

(n + l)(n+2)(n + 3)(n + 4) (1 - (i 2 ) 2 _ 1 , . 

+ 2^ ]? ••••_]' { * 7) 

which, by lemma (ix), must be equivalent to P n /r n + 1 . 

Since (y + V — Is)" is a solution of equation (1) of the n-th degree and 
homogeneous, we have by lemma (ix) as another solution this expression divided 
by r 2n+1 , 



(y + V — lz) n _ (ff + V — MT i 2ra + l2/ 2 + g a 
r 3»+i — a; 2,l+1 L 2 a; 2 

, 2« + 12n + 3 1 (^ + ^) 2 2n + 1 2n + 3 2n + 5 1 (y 2 + s 2 ) 3 _l "1 ajn 
+ 2" 2 2! » 4 ~ 2 2 2 3! x« f •••']> W 

which may, by putting ^(jc) = l/x 2n+1 , be written 

(2 , + v-^r ,)« [ % (.) - ^±^ ^' (-) + ^2l{^FI^+^ % IV (-) 



2 6 3!(n+ 2 l)(n + 2)(» + 3) ^^ + ' * * ' J* (49) 



Putting ^ (as) = shih a;, we have as a solution 

Binn»(y+v l a; |^i 2 2 (n + 1) ^ 2* 2! (« + 1)(» + 2) 



(*/ 2 + z a ) 



2\3 



2 6 3!(rc + l)(n + 2)(ra + 3) 



+ 



..]. (50) 



Since (y + V — lz)*= (cos »<£ + "/ — 1 sin n$) (*/ 2 + a 2 )" /2 , this solution may 

be written 

c n sinh x (cos n$ + V — 1 sin «$) J"„ (®), (51) 

where c„ denotes a constant dependent upon n. 
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If 4" (#) = as" -8 , then (49) becomes 



{y + v—i *y .- [i - (n "gg , - 1 ' ) - 1) 



y 8 + s 2 

a 2 



+ 



(w — s)(n — s — l)(n — s — 2)(n — s — 3) (y 8 -f a 2 ) 2 



2 4 2!(s + l)(s + 2) a; 



•••], (52) 
/• i / 7 \s r \m-sI~! ( n — s)(n — s— 1)1 — a 8 

=(y + v - 1 2 ) s (r^)» s |^i-^ — 2 2 ; ( ; + 1) — —^f~ 

( n -s){n-8-l)(n-8-2){n-s-Z){l-^Y 1 

+ 2*2!(s+l)(s + 2) ^ -.-..J. [bi) 

Since y % + z 2 = r 2 — x % , (52) can be broken up into terms containing the same 
powers of the variables as do the terms of (31); therefore their coefficients are 
proportional and so (53) may be written 



c«.. r n (cos s<p + */ — 1 sin s$) T<f>, (54) 

where c„ i8 denotes a constant dependent upon n and s. 
It is obvious that if in (46) divided by r", or 

P M (cos0) = cos»0[l-^^tan*0 + 

we put 6 = wjn, and then cause n to become infinite, we shall have 

S'P ) ,(co S 3 = l-^ + ^-....=/ (4 (56) 

Similarly, if in (53) or 

r n ~ s {y + */~=l a) 8 cos"- 8 6 [l - ^^ ^h^"*" 1} tan 2 

(n— «)(n-a-l)(»— a — 2)(n — s — 3) 4fl "j , . 

+ 2*2!(s+l)(s + 2) tan »-.... J (57) 

we put 6 = ox/s, and then cause * to become infinite in the part which is multi- 
plied by r n - g (cos$ + */ — 1 sin<£>) s , we shall have for this part 

w '\. 1 ~2»(s+l) + 2*2!( S + l)(* + 2) -■■••] = C M«)> ( 58 ) 

where c' g denotes a constant dependent upon s. 
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According to lemma (ix), xfr 3 is a solution; it is, in fact, the velocity- 
potential defining the motion of a sphere in an infinite mass of water, x/r 3 is 
the kind of solution satisfying condition (5). If in 

x - 1 _ IjJ y 2 + z 2 1.3.5 J_ (y 2 + z 2 ) 2 _ , , 

r 3 ~x 2 2 x* "*" 2.2 2! x 6 "" K > 

we put X (*) = A, then jjr" (a) = ^ , X IV (*) = 2 ' - '* ' 5 > • • • • , we have 

which, being of the same form as (42), leads to no solution additional to those 
which might have been obtained from 1/r. 

In studying Bessel's functions, it is more natural to begin with solutions 
of (2) than with solutions of (1). We have as an illustration the physical problem 
of the vibration of a circular drumhead, or of wave-motion in a circular sheet of 
water, the periodicity of the phenomenon suggesting circular time-functions. 

By lemma (viii) we have as a solution of (2) 

1 _1 , l y 2 + g 2 4 , 1.3 (y« + s 2 ) 8 1.3.5(y 2 + z 2 ) 3 

*/x 2 —(y 2 + z 2 )~x~ r 2 x 3 " r 2.4 x 6 " t "2.4.6 x 7 -f . . . . L»>i; 

Putting x (a) = 1/x, this becomes 

z(*) + l^(y z + z z )z"(*)+ 1 ^lj ] (y z + * 2 yz iv (*)+ (62) 

Suppose, now, that x, y, z are replaced by t, x, y in this expression and in 
equation (2). Let 4> (t) = cos t ; then 4'" (0 = — cos t, 4* lv (*) == cos t, . . . , 
and we have as a solution of (2) thus altered 

™at[l-±A** + y Z )+2^(x z + y z y-....]=co S tJ (w), (63) 

where m 2 = x 2 -\- y 2 . 

If -^{x) = x n we have a n P n f-j as a solution of (2), a 2 denoting x 2 — 

(y 2 + 2 2 )- This solution could have been obtained from r n P n (-) by lemma (iv). 
By lemma (viii), 

(64) 



or 



t 2 — {x 2 + y 2 + z 2 ) t 2 — r 2 
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is a solution of (2) when x, y, z, w are replaced by t, x, y, z. In 

1 _ 1 , r 2 r* r 6 , * 

= 75+ 71 +76 + 78+ •••• ( 65 ) 



* 2 — r 2 t* • t* • t® ' t 
put £ (tf) = l/£ 2 and it becomes 

Let •>// (£) = sin £; we then have as a solution of (2) 

V 4 = sin t ™H. (67) 

r 

The solutions obtained in the next section are also of the kind specified in 
lemma (v), as will be shown presently. 

§ 5. Solutions Obtained by an Application of Taylor's Theorem. 

Suppose q>(ax + (3y + yz + ) to be expanded in powers of (iy + yz 

-+- . . . . thus, 

*(««) + 0y + y»+-- $ (a x) + (ll+^+^y ¥ (a x) 

+ (^ + y«+—)V (aa)+ ,... f (68) 

the accent denoting differentiation with respect to x. (If the accent denote 
differentiation with respect to ax, then, of course, a would not appear in the 
denominator.) 

Equation (1) will be satisfied by this expression if 

a 2 + £ 2 + y 2 + ....=0; (69) 

and equation (2) if 

a 2 — (£* + /*+ ....)= 0. (70) 

For, differentiating the coefficient of <p (,,) (ax) twice with respect to each of 
the variables y, z, . . . . and adding the results together, we have 

P z (/3y + yz + ....)-» y z (/3y + yz + ...Q- 8 

(r — 2)!a" "*" (v — 2)!a* " 1 ~ 

The term 

(Py + y Z + ....)- » 
(v — 2)! a' 

twice differentiated with respect to jc becomes 



rn^ ♦""•(a*) 



a'^y + y. + ....)- w (gg)> 

(v — 2)!a" T v y 
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Thus it is seen that if q>{ax + fiy + yz + ) be substituted in (1) or (2) 

the coefficient of the resulting f w (ax) will be zero, and so (1) or (2) will be 
satisfied when the relation (69) or (70) obtains. Or what amounts to the 
same thing, the given (y, z, )-coefficient of <g> ( *- 2) (a x), added to the 

( ~— -2 + o— f + • • • • )- derivative of $ (y) (ax), is zero, and so this comes under 

lemma (v). 

[It may be noted here that (68) is a solution of the equations 

d*V ,d n V ,d n V , _ A ,,__. 

-^ + dr + '^ r + "" °' (71) 

d n V /d n V d n V \ _ 

~dx^ ~~ \dr + a^~ +••••)- °> ( 72 ) 

provided that, in the respective equations, 

a" + /3* + y*+ .... =0, (73) 

a n — (P n + y n + ) = 0. (74) 

This can be shown by differentiating the coefficient of $ M (ax) n times with 
respect to y, z, . . . . , and adding the results together, and by differentiating 
the term 

(£y + yg+ •■•.)-» <,_„) , * 
(v-n)\a"- n * [aX) 
n times with respect to a?.] 

The whole of the expression $ (a* + fty + yz + ) satisfies (1) or (2). 

Consequently, if it can be thrown into parts linearly independent of one another, 
each such part must be a solution. 

For separating the solutions of (l) just referred to, it is convenient to make 
use of a complex variable ax + iby-\-jcz + Jcdw-\- .... of^> dimensions, 
where the original units are 



1, i = V-l, y=V — 1, 7e = V-l, ...., (75) 

and where 

a 2 = b z + c z + d* + (76) 

Let p denote the dimension of the complex or the number of original units ; 
then the total number of units in the (1, i,j, . . . .)-algebra is 

P + c 9 (p~-l) +c 3 (p-l)+ .... + <V,_ 1 (jp-l). (77) 

Here c z (p — 1) means the number of combinations of p — 1 things taken two 
at a time. 
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Similarly, for separating the solutions of (2) it is convenient to make use 
of a complex variable a.x + iby+jcz-\-kdw + .... of p dimensions, where 
the original units are 



1, i = V + l, j=V + l, k = V + l, ...., (78) 

and where 

a 3 = b 2 + c 2 + d* + 

The parts of (68) thus separated may in the first instance be denoted by the 

Italic symbols 

B u %B t , jBj, ijBij, ...., (79) 

and in the second by the Roman symbols 

R 1? iR lf jRj, ijR n , (80) 

If jP = 2, the two solutions of (1) are the ordinary X, T Of the function- 
theory; the two solutions of equation (2) are somewhat analogous functions 
which may be denoted by Roman X, Y. 

In the case p = 3, we have as solutions of (1) 

Bi + iBi +jB j +ijB ij = <p(ax + iby +jez) — $ [ax + i(by — ijcz)] ; (81) 
B, =^(ax)- b ^l+f z \ "(ax) 

+ " 4{a i ~ <l> lv (ax) -...., (82) 

A^^M- ^' + yy'V H 

+ b'y>+10b'c*y^ + 6bc*^ r{ax) _ ^ (g8) 



d cz A/ , x c 3 z 3 + 3cb 2 zy z .... , N 
j = -£*( ax ) 3!^ -f( aa! ) 

+ c>z>+10cWj>+5cV Z y* r{ax) _ f (84) 

Bii = h ^*»{ax)- h * G y*i+*/y*\ " { ax) 

+ h % *lf* + *W* + 1 *** *"{**)- (85) 

.'. B 1 -B tj + ^^^{B i + B 3 ) = ly\r' 2{ ^p^^(ax). (86) 
In all these cases the accent denotes differentiations with respect to x. 
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For the corresponding solutions of equation (2) we have 
K x + i R s + j Rj + i j Rii = 4> («* + i by + j cz) = <p [(ax + i 6y 4- ij cz)]. (87) 

The R's may be obtained from the preceding expressions by changing each 
minUs sign into a plus. 

If v = 00 , 

JR, + R v = cos pJ^p£_« ^] <?> («x), (88) 

R, - R tj = cos p*±« |J ^ (ax), (89) 

^ + i2, = sin [-^- Z ^] <?> (« «), (90) 

R t - J2, = sin [^^ /-] <?> (a4 (91) 

For corresponding combinations of the R's, replace the symbols cos, sin 
by cosh, sinh. 

As already stated, (68) is a solution of (1) or (2) of the form referred to in 
lemma (v), and so the linearly independent parts R lf R it R j} R ijf .... are 
severally solutions of that form. 

A lemma similar to (v) is applicable to solutions of (71) and (72); the 2's 
throughout equation (5) must then be changed to n's. 

§ 6. Applications of Complex Variables of Two or More Dimensions. 

Instances of the use of the ordinary two-dimensional variable x + iy, or 

x — iy, and the variable \ogr-\-iB (where r 2 = x 2 + y 2 , $ = tan -1 " J for 

finding or expressing solutions of equation (1) need not here be given, since 
they frequently occur in cartography, in works on hydrodynamics, electricity 
and magnetism, and other branches of mathematical physics. 

A few examples of such solutions of (l) and (2) will now be given as may 
be expressed by aid of the other complex variables mentioned in the preceding 
section. The examples show that solutions of this kind, or aggregates of them, 
may have physical applications ; also, that they are in many ways analogous to 
solutions made up of the ordinary X- and F-functions. The conception of these 
complex quantities may not be of great importance in the obtaining of solutions 
satisfying prescribed conditions; but when solutions have been once obtained, 
this idea often facilitates their inspection. 
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Example (1). Any function $ of the variable x-\-'\y, where i = V + 1, 
gives two solutions of equation (2) ; that is 

<p(x + iy) = 'R 1 + i~R i (92) 

gives Rj , Rj as solutions. 

If for x we put xt, and for y, x, then R 1? Rj are solutions of 

3 2 V 3 2 V A 3 2 V 2 8 2 V 

For a vibrating string, the displacement V should be periodic with respect to t. 

7t 71 

$ may be put equal to A x cos -j or ^4 X sin -,- , where I denotes the length of the 
cord. Suppose $ = A 1 sin -y ; then 

R 1 + iR, = A 1 sin|(xi + ia;), (94) 

and so 

» „:„"#'„ "a 5 t> a 7t?ci{ • 7ta3 ,__* 

x = A x sin — p- cos -j- , Ri = A x cos — p sm -=- . (95) 

3 R 
The displacement Rj becomes zero when x = 0, a; = ?; and -«— ' = when 

at 

£ = 0. The problem of having the string initially at rest, and conforming to a 
prescribed curve whose equation is 

y or Rj =/(«), 

generally requires for its solution an infinite number of terms like R,. That is, 
when t=0, y is of the form 

A. 7105 . a • 27135 . i • 37105 , ,„„, 

ism-y- + A 2 sin— = |-A 3 sin— j- + (96) 

In this case the function <£> is really 

Ri + iR 4 = Aj sin | (xt + ix) + A 2 sin ^ (xt + ix) 

+ A 3 sm^(xt + ix) + (97) 

Example (2). Using the variable ax + iby + jcz, where 1 = ^ + 1, 
j = s/ + 1 and, of course, a % = 6 3 + e a , any function 

q> (a x + i b y + j oz) = ^ + i Rj + j R 3 + i j R tj (98) 

gives four solutions of (2). If for as we put xtf, for y, x, and for z, y, the R's 
are solutions of the equation 

3 2 V 9 /3 2 V . 3 s V 



52 



(& + &) = •• <••) 
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Giving to 5, c, the values ^— , -j— , and so to a the value .J ^ + ^ n, 
and making 4> = cosine, we have 



-c . lm, 2 , n 2 . mnx . nny 



■r, . m z , n z mn 

Rj = COS**** Nnr + To cos — =— 

\ A? h% K 



x nny 
- cos ,y 



2 



(100) 



h 

Functions like these are suitable for representing the free vibrations of a 

rectangular or square drum-head, or for representing the motion of long 

stationary waves (i. e. } oscillations) in a rectangular or square tank of uniform 

depth. See example (12) below. 

Example (3). Taking as variable ax + iby -\-j cz, where i = V — 1, 

j = V — 1, a? = b z + c z > and letting 

<}>(ax + iby+jcz) = e -<«> + «»+*•* (ioi) 

we have 

R tj = e~ az sin by sin cz. (102) 

If for b, c we write 

(2m + l) „ (2n + l). 



and so for a, 



"V"*' — *r~~* 



I (2m+ l) 3 (2n + l)« 

\ Af + m 



"2 

we have 

p _ _^V^±^+^® . 2m + 1 . 2n+l , 1AO » 

i%=e wx v m "s sm — • — 7t^ sin — =- *z. ( 10S ) 

"i % 

This solution of equation (1) is applicable to the potential at points within 
an infinitely long rectangular prism, where the potential at the sides is zero, 
and that at the base is a given constant.f 

Example (4). Taking x + iy cos/3 +jz sin/3 as the variable, and e~ i() 
as the function, we have 

$(ax + iby+jcz) =fl -«o»+«iroo.0+^ B ta» ( 104 ) 

which may be written 

e ycoal3 - iJ * ainf> co8x — ismx; 

.-. R 1 + ij R tJ = e y«»i>-v>»i> cosx, ( . 

iRt + jRi = — ie y °°*f>-v* **? sin x { ' 

*See Byerly, "Fourier's Series, etc.," pp. 126-134. + Byerly, pp. 138, 139. 
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are solutions of equation (1). If we do not care to distinguish between the 
1-terms and the */-terms, or the i-terms and the /-terms, we may put ij= ± 1 
and j = ± i in the above results. Multiplying by suitable time-factors and 
combining, we obtain as a solution of (1) 

e VOOBli -* eial> cos(x—xt). (106) 

This denotes the velocity-potential for a train of "edge" waves advancing 
in the aj-direction, the shore sloping at an angle /?. 
This solution is due to Stokes. * 

Example (5). If a 2 = 1, 6 2 = -, c 2 = -, and <p = cosine, we have 

<p (ax + iby +jcz) = cos (a; + ^~y + -^ s) ; (107) 

.". R, = cos x cosh -^= cosh — §= . (108) 

1 V2 */2 v ; 

This expression, when multiplied by the time-factor, denotes the velocity-potential 
of a standing wave in a V-shaped canal whose sides make angles of, 45° with the 
vertical. This type was discovered by Kelland.f 

— R,, = cos x sinh -%= cosh — — 
denotes the velocity-potential of an asymmetrical mode, noticed by Greenhill. J 

k/ ^ 1 

Example (6). Considering now two variables x-\- ix, x -f- i — jr- y +j -z, 
whose coefficients evidently satisfy the relation a z = b z + c 2 , the R's of 

</3 



A cos (<e + iz) + Ccos (x + * -x- y + /o z j (1° 9 ) 



and of 



5 sin (a; + iz) + D sin f a; + * — - y +j^ zj ( 110 ) 

are evidently solutions of equation (1); so will be the R t of the former added to 
the Ra of the latter. The result is 



cos a; 



A cosh 2 + B sinh 3 + cosh -— - y ( cosh - z + D sinh - z j . (Ill) 

* Lamb, "Hydrodynamics," 3d Ed. (1895), p. 431; 3d Ed. (1905), p. 434. 
f Ibid., 3d Ed., p. 433; 3d Ed., p. 434. 
t Ibid., 3d Ed., p. 435. 



412 Harkis: On Harmonic Functions. 

This expression, when multiplied by the time-factor, is the velocity-potential for 
a type of standing waves in a canal whose sides are inclined to the vertical at 
an angle of 60°. A, B f C, D are determined by the boundary and surface con- 
ditions. This comprehensive type was discovered by Macdonald.* 
Example (7). If 

$ (ax + iby +jcz) = log (ax + iby +jcz), (112) 

then 

R i -B j = t^Qy^-) (us) 

must be a solution of equation (l), provided, of course, that a z = b z + c z . 
Example (8). If in 

ty (ax -f- iby + Jcz) = (ax -\- iby + jcz) z (114) 

we make a = v'2, 6 = 1, c = l, then 

^ = 2^-^-^, B i =2V~2xy, 1 

i^ = 2V2a:z, B tJ =2yz. \ ( ' 

If we make c = 0, and so a = 1, 6 = 1, then 

B 1 = x z — y z , Ri = 2xy. (116) 

These are rational integral functions of degree 2 satisfying equation (1), and so 
are harmonic functions. 

Keplacing the exponent 2 by 3, 4, .... , rational functions of those degrees 
are obtained which satisfy equation (1). This fact has been noted by Laurent, f 

For the exponent 3 we have 

Ri + iBi +jBj + ijBv = V^(2x 3 — Zxy z -Zxz z ) 

+ i(6x z y—3yz z — y 8 )+j(6x 2 z — 3y z z — z 8 ) + ij6V~2xyz. (117) 

But 

Bt= V2[2x 8 — 3x(y z + z 2 )], (118) 

having circular symmetry about the z-axis, is, like 2x s — y z — z z , a zonal 
harmonic. 

* See Lamb, "Hydrodynamics," 2d Ed., p. 484; 8d Ed., p. 426. 

t "II est bon d' observer que, si 

a 2 + /3 2 + 1 = 0, (ax + fiy + «)» 

satisfait a A 2 u = et (a sin 6 cos iji + (3 sin sin ip + cos 6) m est une fonction spherique du degre' m. Une 
somme de 2 m + 1 semblables fonctions sera la fonction spherique la plus generate. " H. Laurent, Traiti 
d' Analyse, Tome VI, pp. 297, 298. 
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If we raise V 2~a: + iy +jz to the fourth power, the harmonic function B x 
thus obtained no longer possesses circular symmetry about the a;-axis, although 
symmetrical in y % and z z ; for then 

B 1 = 4x i — 12x*y i —12x*z 2 + Qy*z % + y i + z i . (119) 

Suppose we raise the variable y + iz to the fourth power; its B t or B t 
must satisfy (1). [Cf. lemma (iii)]. 

Adding its B lf = — 6y 2 z z + y* + z*, to twice the preceding B lf we have 
a zonal harmonic of the fourth order 

8x i —24:x z (y 2 + z 2 ) + 3(y 2 + z*y. (120) 

Were it desirable to do so, it seems probable that the zonal harmonic of any 
order could be obtained by combining the i2/s of this power of variables of two 
and three dimensions. 

Example (9). Still considering the case where p = 3, it will be noticed 
that the complex quantity ax + iby -\-jcz contains one effective variable 
parameter. Let a = 1 ; then by (76), the quantity may be written x-\-iby 
+ jVl — b % z. Hence, whenever in the expansion of $ (a? -f i 'b y + j V 1 — b z z) 
it is possible to pick out the entire coefficient of any power of b, that coefficient 
is a solution of (1). This can be easily done for 

(x + iby+jWl—b z z) n , (121) 

where n is a positive integer. The solutions thus obtained are 2n + 1 in number, 
as may be seen upon considering the four kinds (82)-(85), the first containing 
no odd powers of y or z: 

From B lt 
" 7? 



Rj i 



B t 



j> 



n odd 


n even 


i + ==i 


l +5 


n + 1 


n 


2 


2 


i + V 


1 + .-. 


n—1 


n 


2 


2 



(122) 



Total, 2rc+l 2n+l 

The coefficients of the various powers of b (including those of b°) in the 
several kinds of terms may therefore be taken as the 2 n + 1 independent 
solutions of (1) of the n-th degree. 
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Now since the solution r* P n is an homogeneous polynomial of degree n 
where y and z enter by even powers only, and which satisfies (1), it follows that 
it can be made up of linear combinations of the solutions of the first kind 
obtained by the expansion of (121). Hence, the surmise made at the end of 
example (8) turns out to be correct. 

We can assume that if any solution V be a regular function of x, y, z within 
a finite domain surrounding as', y', z', it can be expanded in powers of X, Y, Z, 
where x = x 1 + X, etc., and the series will there be absolutely and uniformly 
convergent. The various powers of the variables x, y, z will give rise to homo- 
geneous polynomials of all degrees (not exceeding n) in X, T, Z. Bach is a 
solution of (1) because, by hypothesis, the whole series is a solution and because 
of the homogeneous character of (1). From what has just been shown, it fol- 
lows that for any such polynomial of degree n there are at most 2 n + 1 inde- 
pendent solutions. Hence, if the expansion be carried to the n-th powers of 

X. Y, Z, we obtain at most 2 2?> -f- 1 or (n + l) 2 solutions.* 

Example (10). Taken in connection with the foregoing section, it seems 

highly probable that - can be expressed in a similar manner by using the B x of 

an infinite number of terms of the form 



ax-{-iby+jcz 



(123) 



We have 



It? f 1 i- 1 I V2 1 

2 1 L2(x + iy)~ r 2(x+jz)^ V2x + iy+jzA' 

_1T 1 l y 2 + 3 2 ,1.3 (y* + z*Y 1.3.5(y* + z*) 3 1 

~ xl 2 x 2 1 "2.4 x i i "2.4.6 x 6 T ''J' 



(124) 



which agrees with the value of 1/r (41) to the number of terms here given. 

It may be noted that a term of the form (123) becomes infinite only when 

the variable becomes zero or lies upon one of the two lines in the yz-plane 

whose equations are 

x = 0, by=±cz. (125) 

* Taken mainly from an important paper (q. v.) by Professor E. T. Whittaker, entitled "On the Partial 
Differential Equations of Mathematical Physics," MathematUche Annalen, Vol. LVII. He appears to have 
been the first to show that (131) gives the general homogeneous polynomial solution of (1). 
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For, (123) when rationalized in denominator gives 

ax — iby — jez 
a 2 x 2 + b 2 y % + c 2 g 2 — lijbcyz 

__ (ax — iby — jez)(a 2 x z + b 2 y 2 + c 2 g 2 + 2ijbcyz) (i 9(i \ 

(a 2 x 2 + b 2 y z + c 2 z 2 ) 2 — 4b 2 c 2 y 2 z 2 ' ' ( ' 

The denominator of this expression vanishes whenever 

(ax) 2 + (by=pcz) 2 = 0. 

Example (11). Consider the i^ terms of (123) or (126). They are 

ax(a 2 x 2 -\- b 2 y 2 + c 2 z 2 ) /i97^ 

(a 2 x 2 + b 2 y 2 + c 2 z 2 ) 2 — ±b 2 c 2 y 2 z 2 ' *• > 

Putting a = 1 and eliminating c by means of the equation — a 2 -J- b 2 -f- c 2 = 0, 
the fraction becomes 

g[a» + 8' + 6'(y»--g')] M9R x 

(a; 2 + a 2 ) 2 +2& 2 (» 2 2/ 2 — a 2 g 2 — # 2 g 2 — a*) + J*(y» + z 8 ) 8 ^ ; 

and may, for brevity, be written 

x\x 2 + z 2 + (y 2 — z 2 )b 2 ] 
p + qb 2 + rb* 

= xlx 2 + z 2 + (y 2 -z 2 )b 2 ][^-^b 2 + 1 ^+^)b^-...^, (129) 

a development in descending powers of (x 2 + g 2 ) 2 . 
The coefficient of b° is 

"FT? > (130) 

the coefficient of b 2 is 

either is a solution of equation (1). 

Example (12). For phenomena periodic in time and involving the differen- 
tial relation 

d 2 2 /3 2 , d 2 



^-'Gk + sW' < 132 ) 



di" \ox" oy 

<£> should consist of a number of sine or cosine terms or of both. Here suppose $ 
to consist of terms of the form 

cos (x at + \bx -f- jcy). (133) 
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(134) 



If the phenomena be simply harmonic in time, then the square of the a of 
each term of the function will be the same. The R's become 

R x = cos xat cos bx cos cy + cos xat cos Vx cosc'y -4- ••••, 
— Rj i= sinxat siu 6a; cos ay -\- sinxat sin 6' a; cos dy + ••••, 
— Rj = sin xa< cos bx sin cy + sin xat cos 6'x sin c' y ■+• . . . ., 
— Rij = cos xat sin 6a; sin cy + cosxatf sin 6' x sin c'y + • • • •, 
where a 2 = b % + c z = b' z + c n = ..... 

As already noted, functions of the form cos xat sin bx sin cy or cos*a£ 
cos bx cos cy are suitable for representing the vibration of a rectangular drum- 




FlGt. 1. 



head or the oscillation of a shallow rectangular body of water. It remains to 
suggest here that by taking a sufficient number of similar terms it would probably 
be possible to thereby represent these phenomena in triangular, circular, or any 
other shaped areas. For instance, the solution given by Lame* for the vibration 
of a membrane whose form is that of an equilateral triangle will be seen to 
consist of terms of the form (134) as soon as his trilinear coordinates are replaced 
by Cartesian coordinates conveniently chosen. We shall here confine ourselves 
to the oscillations of shallow bodies of water. 

One mode of oscillation for a 45° triangle is shown in Fig. 1. The vertical 
displacement, £, is proportional to 

7txt / 
( c 



cos 



h \ 



nx , nv\ 
cos x + cos£), 



(135) 



'Lesons sur la theorie mathSmatique de Pelastieit.6 des corps solides," p. 132. 
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which is of the form 

cos xat [cos b x cos c y + cos V x cos d y], 

where c 2 = b n = 0, b 2 = c' 2 = a 3 , and x 2 = g X. depth of water. 

Had the coordinate axes been taken along the legs of the triangle, the 
above displacement would have been written as one term, viz., 



which is of the form 



<, nxt nx ny 

2 cos cos — t=— cos -y , 

A V 2 7t V 2 A ' 



2 cos xaf cos bx cos c^, 



(136) 




Fio. 2. 

where b z = c 2 = -a 2 . (Of. Example 5, above.) 
■a 

If the above expressions for vertical displacement be expanded in powers 
of x and y, the results are the same to two terms, viz., 

2 cos xat (l — a 2 — j" 2 y * ^\ . 



(137) 



Considering now a square composed of four triangles, arranged as in the 
figure, we see that for this mode of oscillation there is a single nodal line con- 
sisting of four straight lines in the form of a square. Upon differentiating 
normally to any boundary of any of these triangles, or to the lines bisecting 

their right angles, it is easily seen that ^- = 0. 
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One mode of oscillation for an equilateral triangle is shown in Fig. 2. 
Here £ is proportional to 

s/Znx ny -»."i * 



TtxtYn \/3nx 7ty , Ttyl : 
cos T [2 cos -^ cos ¥ | + cos -/J , 

which is of the form 

cos xat [2 cos bx cos cy + cos 6'a cos cf y], 



(138) 



w 



here b 2 = -a 2 , c 3 = - a 2 , 6' 3 = 0, c' 2 = a 3 . (Cf. Example 6, above.) 




Fig. 3. 



This, if expanded, gives for the first three terms 

3 cos xatyi — a 3 — J-*- + « 4 V 2 ' 4' J ' 



(139) 



Considering now a hexagon composed of equilateral triangles as in the 
figure, we see that for this mode of oscillation there is a single nodal line 
resembling the circumference of a circle much more nearly than did the nodal 
line in the preceding case. Upon differentiating normally to any boundary of 
any of these triangles or to lines bisecting their angles at the center of the 

hexagon, it is easily seen that «-^ = 0. 



* Coast and Geodetic Survey Report for 1900, p. 589; also Poincare, "Lecons de Mecanique Celeste," 
Tome III, p. 371. 
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The vertical displacement for the slowest symmetrical mode for a circle 



is proportional to 



cos xat J (a m), 



or 



cos x 



at^i—a — ^2 — + a — p-42 a 2 2 4 2 6 2 + J' ' -* 

That this satisfies the differential relation (132), has been shown under the 
heading "Application of Lemma (v)". 

At the circumference where m = m lf d J (am)/ dm must be zero. Such 
a value of am is am 1 -=- 1.22 n, or a = 1.22 nlza 1 ; that is, ^ being given, 
a becomes fixed, and conversely. If a m = 0.7655 n, J (am) becomes zero 
and so does the displacement £. 

.-. 1.22 — = 0.7655; and so m — 0.628^ 
m 1 

is the radius of the nodal circle. Since the motion is symmetrical about the 

Pi Y 
center, the vectorial angle 4> is not involved ; consequently, ~-^ = everywhere. 

Here the "speed" (i.e., 2 n — period) is xa or 1.22 x7t/m 1} while the 
"speed" for a rectangular trough whose length is m x is xn/m lt x denoting, 
as usual, the velocity of advance dependent upon the depth. The period for 
the circle is, therefore, 1/1.22 = 0.820 times the period for a rectangular 
trough m x in length. 

Certain ratios for square, hexagonal and circular areas are given here : 





Height (?) 

at center -s- 

Average height 

along perimeter 


Area enclosed 
by nodal line -*- 
Area of polygon 


Virtual 
length (A) -4- 
Apothem of 

polygon 


Square (Fig. 1) 
Hexagon (Fig. 2) 
Circle (Fig. 8) 


- 8.000 

- 2.418 

- 2.488 


0.500 
0.899 
0.394 


1.000 
0.866 
0.820 



Comparing expressions (135), (139) and (140), we are led to believe that 
the B x of a sufficient number of terms of the form cos (xat + ibx + jcy) 
could be made to represent this oscillation of a circular sheet of water, just 
as in the preceding example we were led to believe that the solution, 1/r, of 
Laplace's equation could be regarded as the R t of an infinite number of terms 
of the form 



ax + iby + jcz' 
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Moreover, it has been shown that any zonal harmonic of order n can be regarded 
as the .ft, of a suitable number of terms of the form 

{ax + iby + jcz) n . 

§ 7. Differential Relations Connecting the R- Functions. 

If d(B 1 + iB i +jB j + ijB ii ) (141) 

d{ax + iby +jcz) 

has a unique value at a given point, it must be free from the variable quantities 
dx, dy, dz. Whenever the numerator contains the denominator as a factor, 
we have the relations 

lfdB 1 dB i .dB J ..dB ij \ 

l/dR l ,.dR i ,.dR J , ..d B tj \ ,,, A „. 

= ^- ( -5—1 + I -5— 4 +J -jr-1 + %J -5-^ J V (142) 

ib\oy dy — - dy * dy J 

Jc~Vd^ +t ~dz~ +J ~d¥ + iJ ~dJ~J- 

• dRl — dR i = dB J 
adx bdy cdz' 

dB x „ _dRi~ dB u 
bdy adx cdz ' 

dR x _. dRj__ dR l} 
cdz adx bdy' 

_ <L?i = _ Ml = dR ii 
cdz b^y adx 

are therefore the nine conditions for a monogenic function of the three-dimen- 
sional variable ax + iby -\-jcz. 

When the Roman R's and i, j are used, all signs are positive in these 
relations. 

If the variable have more than three dimensions, the additional relations 
can be readily determined. 

For two-demensional variables we have, of course, 

dX dY dX dY 



(143) 



dx dy ' dy dx ' 

ax__aT ax _3y 

dx dy ' dy dx ' 



■ (144) 



